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SECTION  1 


INTRODUCTION 


The  t  '  e  Control  of  Space  Structures  (ACOSS)  Eleven  contract  includes 
three  'Ustitiu  areas  of  endeavor:  Simulations  Exf^nsion,  HALO  Optics,  and 
ACOSu.  This  report  covers  the  work  performed  in  each  of  these  areas  between 
April  and  September  1981.  A  brief  introduction  to  each  of  the  tasks  covered 
follows. 


1 . !  Simulations  Extension 

In  order  c-  generalize  the  HALO-or iented  portions  of  the  Draper  Integra¬ 
ted  ^inv. lat ions  vDIS)  and  simulate  other  useful  surveillance  concepts,  CSDL 
i.i  tiaten  three  major  activities  during  this  reporting  period: 

•' t 1  )  i  .ial  planning  with  DARPA  to  select  surveillance  systems  concepts 

>  simulate  an-  evaluate. 

(.1)  Establishing  the  capability  to  generate  and  manipulate  synthetic 

scenes  in  house  as  -ell  as  support  DARPA  in  standardizing  key  scene 
data  bases. 

(3)  Extending  the  current  D1S  algorithms  to  handle  a  wider  class  of  in¬ 
put  scene  data  sets. 

The  Gemini  concept,  being  developed  by  Aerojet  Electro  Systems,  is  the 
only  system  presently  identified  for  simulation  and  evaluation.  As  first  or¬ 
der  Gemini  models  become  available,  they  will  be  incorporated  into  the  DIS. 

A  "generic  scene  generation"  capability  is  being  established  as  part  cf 
the  DIS  to  extend  its  range  of  applications.  Initially,  the  simulation  will 
contain  a  Jata  base  that  is  a  geometric  representation  of  an  '40  *  40  5cm  re¬ 
gion  in  Southern  California  and  it  will  have  a  spatial  resolution  of 
'100  m.  A  set  of  four  prescriptions  will  allow  this  data  set  to  represent 
typical  scenes  from  the  Soviet  tundra,  the  Arctic,  the  Middle  East,  and 
Central  Europe  with  a  limited  capability  to  generate  clouds  and  superimpose 
them  on  any  geographic  region. 

Photon  Research  Associates  (PRA),  under  subcontract  to  CSDL,  is  develop¬ 
ing  four  standardized  surveillance  scenes  to  assess  various  surveillance  sys¬ 
tems  concepts.  DARPA  will  gain  a  baseline  to  compare  system  performance 
evaluations  for  common  missions  and  a  pr  -cedure  to  identify  the  scene /system 
characteristics  that  are  performance  drivers  from  this  work. 


HALO  Optics 


Three  contractors,  Eikonix  Corporation,  Itek  Corporation,  and  Hughes 
Aircraft  Corporation  (MAC),  have  been  working  on  the  general  deconvolution 
problem  which  is  classified  as  the  HALO  Optics  task  of  this  contract.  Eikonix 
has  worked  on  the  phase  retrieval  problem  whete  aberrations  are  determined 
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from  the  focal-plane  data;  Itek  has  used  a  wavefront  sensor  to  measure  the 
system  aberrations;  and  HAC  h3s  used  "color"  algorithms  on  the  focal  plane 
data  to  obtain  optical  system  correction. 

During  this  reporting  period,  CSDL  developed  a  computer  simulation  to 
generate  noisy  and  aberrated  *>cal-plane  point-spread-function  (PSF)  d*  a 
under  guidelines  from  Eikonix. 

CSDL  a  so  prepared  a  phase  retrieval  test  based  on  aberrations  derived 
from  cryogenic  deformations  of  an  Itek  HALO  mirror.  The  aberrated  PSF  was 
sampled  by  an  8  x  8  array  of  square  detector  elements,  each  of  which  had  a 
full  width  of  2.13  AF.  Three  test  cases  were  run.  Case  1  had  Gaussian  random 
noise  which  had  a  uniform  2-percent  standard  deviation  of  the  peak 
diffraction-limited  signal.  In  Case  2,  the  PSF  was  decentered  such  that  the 
peak  irradiance  of  the  diffraction-limited  point  spread  function  was  moved  to 
correspond  to  a  line-of-sight  error  of  (0.984  A/D,  -0.984  A/D),  where  D  is  the 
diameter  of  the  system  pupil.  As  before,  Case  2  had  2-percent  noise.  In  Case 
3,  there  was  no  line-of-sight  error,  but  the  noise  was  increased  to  5  percent. 

HAC  reviewed  their  work  on  OYSTER,  color  algorithms,  image  moments, 
CORRWAVE ,  and  phase  retrieval  at  DARPA  in  May  of  1981.  CSDL  wilx  prepare  test 
cases  for  the  HAC  color  algorithms  during  the  third  year  of  the  ACOSS  11  pro¬ 
gram. 


1 . 3  Active  Control  of  Space  Structures  (ACOSS) 


Volume  2  of  this  report  describes  the  work  CSDL  has  done  to  investigate 
spacecraft  control  theory.  Each  of  the  six  sections  devoted  to  ACOSS  reports 
on  a  different  aspect  of  that  work. 

Section  4,  "Compensated  Truncation  of  Modal  Models  for  Design  of  Control 
Systems,"  describes  the  selection  process  necessary  in  large  space  structure 
(LSS)  control-system  design  using  a  truncated  finite-element  model.  The  trun¬ 
cated  model  must  be  selected  properly  and  compensated  explicitly  for  control 
and  observation  spillover,  so  the  control  system  designed  through  thi3  mtthod 
can  perform  sacisfactorily  when  implemented  on  the  structure.  Proper  selec¬ 
tion  requires  correct  classification  of  structural  modes  into  "primary"  and 
"oecondary”  modes.  Explicit  compensation  for  truncation  incL.  ies:  placement 
of  actuators  and  sensors,  synthesis  of  the  actuator  and  sensor  influences  once 
they  are  placed  on  the  structure,  and  filtering  of  the  actuator  inputs  and 
sensor  outputs. 

Section  5,  "Ensuring  Full-Order  Closad-Loop  Stability  in  the  Reduced- 
Order  T  da  m  of  Output  Feedback  Controllers,"  builds  on  the  studies  performed 
during  ACOSS  6  that  established  various  conditions  necessary  to  ensure  full- 
order  closed-loop  asymptotic  stability  and  robustness  with  raduced-order  de¬ 
sign  of  velocity  and  displacement  output  feedback  controllers.  Currently,  the 
work  In  tnis  area  concentrates  on  how  to  apply  such  results  to  large  flexible 
space  structures  and  how  to  develop  a  reduced-order  design  technique  that  will 
ansure  full-order  closed-loop  asymptotic  stability. 


The  study  includes  preliminary  development  or  computer-aided  design 
sortware  and  acceleration  output  teedback  control. 


Section  6,  "Design  Freedom  and  the  Implementation  o:  Suboptimal  Output 
Feedback  Control,"  discusses  the  freedom  inherent  in  design.  The  section 
states  that  often  this  treedom  is  sacrificed  purposely  when  simplifying  as¬ 
sumptions  are  made  to  avoid  theoretical  or  computational  difficulties.  Since 
it  is  Qirncuit  to  consider  this  topic  without  referring  to  specific  applica¬ 
tions,  the  section  uses  controller  design  as  an  example  where  work  is  being 
done  to  discover  and  exploit  the  freedom  of  choice  in  design.  Then,  the  sec¬ 
tion  uses  suboptimal  output  feedback  control  as  a  case  study  which  is  relevant 
to  ACOSS  development. 

Section  7.  "Stochastic  Output  Feedback  Compensators  tor  Distributed 
Parameter  Structural  Models,"  presents  recent  progress  on  the  stochastic 
output  feedback  design  problem  tor  distributed  parameter  plants.  The  results 
presented  are  an  extension  of  work  done  under  the  previous  contract. 


The  concepts  developed  are  general  enough  to  apply  to  a  wide  variety  of 
tixed-rorm  compensator  design  problems,  and  current  studies  are  aimed  at 
specializing  the  results  to  the  optimal  output  feedback  compensator  design 
proDlem.  The  procedure  developed  will  be  applied  to  the  design  of  velocity 
feedback  controllers  for  a  vibrating  string.  The  results  or  this  simple  test 
should  provide  insight  into  the  impact  of  various  modeling  assumptions  on  the 
convergence  of  the  design  procedure  described. 

Section  b,  “Large-Angle  Spacecraft  Slewing  Maneuvers,"  further  develops 
work  that  was  reported  in  the  previous  ACOSS  contract.  Specifically,  the  sec¬ 
tion  presents  techniques  for  improving  tne  optimal  torque  profiles  by  allowing 
the  solution  process  to  determine  the  optimal  terminal  boundary  conditions  and 
by  developing  a  control-rate  penalty  technique  for  producing  smooth  control 
profiles.  Several  example  me 'euvers  are  provided  to  demonstrate  the  practical 
application  and  utility  of  the  techniques  presented. 

Section  9,  "Order  Reduction  by  Identification— -Some  Analytical  Re¬ 
sults,”  attempts  to  characterize  control  designs  that  will  guarantee  stability 
using  a  reduced-order  model.  This  kind  of  design  compromise  practiced 
regularly,  but  no  one  has  verified  the  validity  or  such  an  approach. 

The  least  squares  (LS)  method  is  used  in  this  analysis  because  it  is  a 
relatively  robust  identitication  scheme  and  analytical  expressions  for  order 
reduction  already  exist  tor  it.  The  results  of  the  analysis  show  that  a  re¬ 
duced  order  controller  can  be  built  using  the  LS  method  of  identification.  It 
is  planned  to  demonstrate  the  practicality  of  this  approach  on  Draper  Model  #2 
in  the  nea.  future. 
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SECTION  2 


SIMULATIONS  EXTENSION 


2 . 1  Introduction 

The  basic  objective  of  the  Simulations  Extension  Project  is  to  general¬ 
ize  the  HALO-oriented  portions  of  the  Draper  Integrated  Simulations  (DIS)  and 
to  simulate  other  useful  surveillance  concepts.  Three  major  activities  have 
been  initiated  during  the  current  reporting  period: 

(1)  Initial  planning  with  DARPA  to  select  candidate  surveillance 
systems  concepts  to  simulate  and  evaluate. 

(2)  Establishing  an  in-house  capability  for  generating  and 
manipulating  synthetic  scenes  as  well  as  supporting  DARPA  in 
standardizing  key  scene  data  bases. 

(3)  Extending  the  currently  implemented  algorithms  in  the  DIS  to 
effectively  and  efficiently  handle  a  wider  class  of  input  scene 
data  sets. 

Presently,  the  only  system  that  has  been  identified  clearly  for  simula¬ 
tion  and  evaluation  is  the  Gemini  concept  under  development  by  Aerojet  Electro 
Systems.  As  this  concept  is  defined  further  and  first-order  models  become 
available,  they  will  be  incorporated  into  the  DIS. 

The  bulk  of  the  work  performed  on  the  Simulations  Extensions  Project 
during  the  current  reporting  period  has  been  in  the  other  two  areas  of  en¬ 
deavor,  and  this  is  discussed  in  detail  in  the  sections  that  follow. 

2.2  Generic  Scene  Simulation 

The  DIS  is  a  sophisticated  analysis  tool  for  overall  evaluation  and  per¬ 
formance  assessment  of  space-based  surveillance  systems;  it  models  the  mechan¬ 
ical,  optical,  control,  signal  collection,  and  signal  processing  subsystems  in 
detail  in  a  highly  interactive  fashion.  Currently,  the  range  of  DIS  applica¬ 
tions  is  limited  because  suitable  data  bases  appropriate  for  the  problem  to  be 
studied  are  unavailable.  As  a  first  step  in  overcoming  this  limitation  and 
enhancing  the  simulation's  capability  to  respond  to  DARPA  analysis  needs,  a 
generic  scene  generation  capability  is  being  established  as  part  of  the  DIS. 
Initially,  this  capability  will  enable  the  DIS  to  generate  and  manipulate  a 
limited  number  of  synthetic  terrestrial  scene  data  sets  as  a  function  of  raajot 
surveillance  system  and  mission  parameters.  Ultimately,  it  is  planned  to  in¬ 
terface  the  generic  scene  simulation  (GSS)  with  the  Defense  Mapping  Agency 
(DMA)  data  base  and  to  use  this  daca  base  as  the  source  for  scene  data  input 
to  the  simulation.  Draper  has  placed  Photon  Research  Associates  (PRA)  under 
subcontract  to  help  achieve  these  objectives.  Seme  of  the  GSS  features  are 
discussed  in  the  following  paragraphs. 
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The  simulation  user  will  be  able  to  vary  independently  the  position  of 
the  observer  (i.e.,  altitude  and  location),  position  of  the  scene,  position  of 
the  sun  (including  night  conditions),  material  types,  width  and  location  of 
the  spectral  interval  (within  the  limits  of  2.5  and  13  pm),  and  observer 
f ield-of-view  and  resolution  with  respect  to  the  scene. 

Initially,  the  simulation  will  contain  a  data  base  that  is  the  geometric 
representation  for  one  generic  scene.  When  related  to  the  earth's  surface, 
the  geometric  representation  will  correspond  to  an  approximately  40  x  40  kilo¬ 
meter  geographic  region  in  Southern  California.  The  scene  spatial  resolution 
will  be  approximately  100  meters.  A  set  of  four  prescriptions  will  be  pro¬ 
vided  to  enable  this  data  set  to  be  transformed  into  geometric  representations 
typical  of  the  following  four  geographical  regions:  Soviet  tundra,  Arctic, 
Middle  East,  and  Central  Europe.  In  addition,  there  wili  be  a  limited  capa¬ 
bility  for  generating  clouds  and  superimposing  them  on  any  of  the  geographical 
regions. 

The  simulation  will  include  a  data  base  containing  approximately  12  com¬ 
monly  found  terrestrial  materials  as  required  for  assigning  characteristics  to 
each  facet  of  any  of  the  five  geometric  representations  cited  above.  The 
simulation  will  also  include  the  LCWTRAN  atmospheric  model,  a  heat  transfer 
module,  and  a  solar  ephemeris  module  for  modeling  the  effects  of  solar  and  en¬ 
vironmental  heating  on  the  scene. 

2 . 3  Use  of  Standard  Scenes 

During  the  current  reporting  period,  DARPA  began  an  effort  to  standard¬ 
ise  the  scene  data  bases  in  use  by  the  surveillance  community  to  assess  the 
performance  of  various  surveillance  systems  concepts.  The  initial  effort  is 
directed  toward  the  space-based  air-vehicle  detection  (AVD)  problem,  and  the 
baseline  scenes  to  be  used  are  summarized  in  Table  2-1.  Each  modeled  scene 
will  be  an  extension  of  measured  data.  In  addition,  each  of  the  scenes  will 
be  modeled  at  two  viewing  angles,  two  times  a  day,  and  in  four  wavelength 
bands.  PRA  will  perform  the  scene  modeling  effort  for  DARPA. 

Use  of  the  standard  IR  scene  data  base  with  the  DIS  will  give  DARPA  a 
number  of  benefits,  the  most  important  of  which  are  the  following: 

(1)  A  baseline  for  comparative  system  performance  evaluations  for 
common  missions. 

(2)  A  procedure  for  identifying  the  scene/system  characteristics  that 
are  the  key  system-performance  drivers. 

The  following  steps  outline  a  candidate  approach  for  conducting  compara¬ 
tive  systems-performance  evaluations. 

(1)  Identify  a  set  of  standard  scenes  appropriate  for  the  systems  to 
be  compared/evaluated. 

(2)  Divide  the  "full"  set  of  standard  scenes  into  two  subsets: 


Table  2-1.  Baseline  standard  scenes  for  the  AVD  mission. 


Scene  Type 


Multi-Layer  Clouds 
over  Ocean 


North  Canada 
Melt  Lakes 


Snow  Covered 
Alaska  Mountains 


Arctic  Sea  Ice 


Location 


Description 


North  Atlantic 


Low  aiticude  linear  structured 
stratus,  medium  altitude  cumulo¬ 
nimbus,  and  high  altitude  semi¬ 
transparent  cirrus. 


Canada  Coastline 
below  Beaufort 


Low  relief  tundra  with  melt 
lakes  surrounded  by  marsh  land. 
Summer  season. 


Brooks  Range, 
Alaska 


Moderate  relief  mountains  on 
north  coast  of  Alaska  south  of 
Beaufort  Sea.  Spring  season. 


Beaufort  Sea 


Snow  covered  sea  ice  of  varied 
thicknesses  with  cracks,  including 
open  water.  Spring  season. 


(a)  A  generally  available  "public"  set  with  targets  with 
characteristics  and  trajectories  that  are  specified  openly. 

(b)  A  "private"  set  containing  targets  with  characteristics  and 
trajectories  that  are  known  by  only  a  limited  group. 

(3)  CSDL  and  each  systems  contractor  will  use  the  "public"  set  to 
validate  mutually  one  another's  overall  performance/simulation. 

(4)  Once  agreement  is  reached  with  the  "public"  set,  CSDL  uses  the 
validated  simulation  to  generate  a  focal-plane  output  for  scenes 
from  the  "private"  set. 

(5)  CSDL  and  each  systems  contractor  then  independently  simulate  the 
signal-processor  performance  in  target  deCect ion/acquis  it  ion. 

Figure  2-1  schematically  illustrates  this  process. 

2 . 4  Enhanced  Scene  Processing 

2.4.1  Introduction 


As  documented  previousl)  [2-1],  the  present  version  of  DIS  performs,  in 
the  spatial  domain,  the  convolution  of  the  scene  radiance  map  with  the  optical 
system  point-spread  function  (PSF)  to  obtain  (in  the  spatial  domain)  the 
focal-plane  irradiance.  Precomputed  image-plane  convolution  tables  increase 
the  efficiency  of  the  algorithms  for  this  technique.  This  approach  was 


Figure  2-1.  Roedmap  for  comparative  systems  performance  evaluations 
using  standard  scenes  and  the  DIS. 


developed  and  implemented  with  the  needs  of  Mini-HALO  in  mind.  Typical  Mini- 
HALO  missions  can  be  analyzed  and  simulated  using  this  technique  with  reason¬ 
ably  modest  expenditures  of  computer  time. 

In  general,  the  computational  effort  required  to  process  a  given  scene 
using  this  technique  is  a  strong  function  of  a  number  of  parameters,  the  most 
important  of  which  are  as  follows:  Scene  extent  (i.e.,  field-of-view  size); 
scene  resolution;  and  line-of-sight  perturbation  characteristics  (i.e.,  line- 
of-sight  drift  ana  jitter).  Thus,  compared  to  Mini-HALO,  the  surveillance- 
system  simulations  with  large  f ields-of-view  and/or  high-resolution  sensors 
and/or  large  line-of-sight  perturbations  can  require  prohib; t ively  long  com¬ 
puter  time.  The  Gemini  system  discussed  in  Section  2.1,  which  requires  high- 
resolution  scenes,  is  a  case  in  point.  There  are  several  approaches  to  gener¬ 
alizing  the  scene  processing  of  the  DIS. 

2.4.2  Fourier  Trans 

One  approach  employs  Fourier  Transform  techniques.  The  scene  is  ini¬ 
tially  transformed  into  the  spatial  frequency  domain  using  Fast  Fourier  Trans¬ 
form  (FFT)  algorithms.  The  result  of  this  transformation  is  multiplied  tnen 
by  the  .ransfer  functions  associated  with  the  optical  system  and  the  sampling 
aperture  of  the  focal-plane  pixels.  The  steps  up  to  this  point  need  to  be 
performed  only  once.  Then,  for  each  sequential  focal-plane  integration  per¬ 
iod,  the  phase  coefficients  of  the  transformed  scene  (in  the  spatial  frequency 
domain)  are  perturbed  to  account  for  the  effects  of  line-of-sight  jitter  and 
drift.  Then  the  amount  of  data  in  this  perturbed  two-dimensional  scene  spec¬ 
trum  is  compressed  by  a  factor  corresponding  to  the  background-to-detector 
oversampling  ratio  and  transformed  back  to  the  spatial  domain  to  obtain  the 
scene  irradiance  per  detector  at  the  focal  plane  during  the  integration  per¬ 
iod.  For  each  sequential  focal-plane  integration  period,  the  phase  coeffi¬ 
cients  are  updated  to  account  for  changes  in  the  drift  and/or  jitter,  and  the 
process  is  repeated. 


2.4.3  Image  Plane  Interpolation  Approar*- 


Another  approach  to  generalizing  the  DIS  scene-processing  capability 
proceeds  along  the  lines  of  the  present  spatial  domain  scene/PSF  convolution 
approach,  except  the  convolutions  are  performed  selectively.  Basically,  the 
concept  is  to  interpolate  previously  computed  sets  of  focal-plane  irradiance 
values  to  obtain  those  corresponding  to  new  positions  of  the  line-of-sight. 


Indeed,  the  present  approach  requires  that  a  new  convolution  be  per¬ 
formed  at  each  instant  of  time  for  which  the  line-of-sight  changes,  e.g.,  for 
each  new  line-of-sight  jitter  step.  However,  consider  a  scene  (a  subset  of 
which  is  the  instantaneous  field  of  view  (IFOV)  of  the  sensor),  and  suppose  £ 
square  reference  grid  is  superimposed  on  the  scene.  For  convenience,  assume 
the  dimensions  of  a  particular  element  in  this  grid  correspond  to  the  projec¬ 
tion  of  the  sampling  aperture  of  a  focal-plane  pixel.  Further,  suppose  that 
this  two-dimensional  scene  is  ov?rsampled  (with  respect  to  a  focal-plane 
pixel)  by  a  factor  of  M  in  each  direction.  Then  an  element  in  this  reference 
grid  will  contain  M  *  M  scene  elements. 


At  the  start  of  the  simulation  run,  select  the  element  of  the  reference 
grid  that  contains  the  line  of  sight  of  the  IFOV.  Then,  that  line  of  sight 
also  will  lie  somewhere  in  a  square  grid  formed  by  four  scene  elements.  If 
the  four  sets  of  focal-plane  irradiance  values  are  computed  that  correspond  to 
the  line  of  sight  coinciding  with  each  of  the  four  scene  elements,  then  the 
focal  plane  irradiance  values  for  the  actual  line-of-sight  position  may  be 
computed  by  two-dimensional  interpolation  in  these  four  sets  of  irradiance 
values. 

Consider  the  line-of-sight  position  at  the  next  time  step,  and  assume  it 
was  displaced  by  drift  or  jitter.  Then,  if  the  new  line-of-sight  position 
lies  within  the  same  square  scene-element  grid,  the  new  set  of  focal-plane  ir¬ 
radiance  values  may  be  obtained  by  interpolation  into  the  previously  computed 
four  sets  of  irradiance  values,  and  no  new  convolutions  are  necessary.  If  the 
new  LOS  position  lies  within  an  adjacent  square  scene-element  grid,  then  two 
(or  at  most  three)  new  convolutions  would  need  to  be  performed  before  the 
interpolation  could  be  conducted.  However,  since  the  end  result  required  is 
the  net  focal-plane  irradiance  per  pixel,  at  most  M  *  M  convolutions  would  be 
required  to  handle  all  LOS  perturbations.  By  selecting  the  interpolation 
scheme  appropriately  (e.g.,  bilinear  or  bicuoic  spline),  it  may  be  possible  to 
use  a  coarser  grid  of  background  samples,  such  as  every  other  one,  and  hence 

2 

reduce  the  total  number  of  convolutions  correspondingly,  possibly  to  M  /A.  At 
each  step  in  the  process,  the  involutions  computed  for  the  scene-element  grid 
points  are  saved  and  aie  aval  le  for  use  in  a  later  run.  Thus,  in  a  set  of 

runs  involving  the  stitne  scent,  the  convolutions  would  need  to  be  performed 
only  for  the  first  run. 
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SECTION  3 


HALO  OPTICS 


3.1  Introduction 

In  the  current  CSUL  HALO  program,  the  general  deconvolution  problem  is 
classified  as  the  HALO  Optics  task.  Three  contractors  nave  been  working  on 
tins  problem.  Eikonix  Corporation  has  worked  on  the  phase  retrieval  problem 
where  aberrations  are  determined  trom  the  local-plane  uata.  Itek  Corporation 
uses  a  wavefront  sensor  to  measure  the  system  aberrations.  These  errors  are 
decomposed  into  mirror-figure  errors  whicn  are  then  corrected  with  actuators 
on  the  mirrors.  Hughes  Aircrart  Corporation  (HAC)  nas  worked  on  the  problem 
ot  opticai-system  correction  using  ' color"  algorithms  on  the  tocal-piane  data. 

CSUL  supports  DAKPA  on  HALO  Optics  technology  oy  preparing  tests  lor  tne 
deconvolution  problem,  evaluating  the  results,  and  assessing  the  performance. 
To  date,  only  software  tests  have  been  prepared  tor  Eikonix,  although  a  hard¬ 
ware  test  is  suggested  in  this  report.  Eor  Itek  and  HAC,  only  software  tests 
are  pianned  in  the  current  program. 

In  this  report,  a  recent  test  is  described  that  was  prepared  tor 
Kikomx's  phase  retrieval  algorithm.  The  preparation  of  the  test  and  the  re¬ 
sults  to  date  are  discussed  herein.  It  is  pointed  out  that  the  aberrations 
retrieved  by  Eikonix  do  not  represent  a  good  estimate  ot  the  actual  aberra¬ 
tions.  When  the  retrieved  aberrations  are  subtracted  rrom  the  actual  aberra¬ 
tions,  i.e. ,  when  a  correction  is  made,  the  residual  aberrations  are  much 
worse.  This  tact  nas  been  reported  to  Eikonix.  Tney  have  been  given  the  ac¬ 
tual  aberrations  to  see  if  they  can  determine  any  error  and  it  they  can 
improve  their  system  performance. 

CSUL  nas  applied  this  test  to  the  image-sharpening  algorithm  with  very 
encouraging  results.  The  Strehi  ratio  ot  the  image  increased  from  U.2o  to 
U.b2  with  one  iteration  oi  correction. 

A  nardware  simulation  should  be  set  up  to  determine  the  limits  or  ap¬ 
plicability  oi  tiie  phase-retrieval  algorithms.  A  procedure  by  which  this  may 
be  accomplished  is  outlined  in  this  report. 

Comments  on  a  briefing  to  tlAC  also  are  reported  herein. 

3.2  Phase  Retrieval  Test 
3.2.1  Phase  Retrieval  Technique 

In  the  Eikonix  phase  retrieval  algorithm,  aberrations  are  estimated  iron 
the  focal  plane  data  by  an  iterative  process  in  wmch  a  merit  function 
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Q,  -  /  j  t(p)  -  x(p)  2  dp 


or 


Q2  -  /[l(r)  -  I(r)  ]2  dr 


is  minimized.  Here,  x(p)  and  t(p)  are  the  aberrated  and  trial  optical  trans- 
rer  tunctions  (OTF),  respectively,  and  I(r)  and  l(r)  are  the  corresponding 
point-spread  tunctions  (PSF).  The  function  I(r)  is  the  measured  focal  p_.ine 
distribution.  I(r)  and  r(p)  are  related  by  a  Fourier  transform. 

3.2.2  Test  Preparation 

3. 2. 2.1  Shortcomings  of  Previous  Tests 

In  past  tests,  the  PSF  data  was  sampled  by  point  detectors  with  a  spac¬ 
ing  ot  XF/5.  Here  X  is  the  wavelength  of  the  object  radiation  and  F  is  the 
tocal  ratio  (ttf)  ot  the  HALO  Optical  system.  Also,  any  noise  in  these  sampled 
data  was  not  considered.  Thus,  the  finite  size  ot  the  detector  elements  and 
system  noise  were  not  included  because  Eikonix  algorithms  had  not  matured  to 
handle  these  aspects. 

3 . 2 . 2 . 2  New  Test 

As  the  Eikonix  algorithm  developed  turther,  CSDL  improved  its  software 
to  overcome  the  shortcomings  ot  the  previous  tests.  Under  guidelines  from 
Eikonix  on  the  size  ot  the  detector  element  and  array,  CSDL  developed  a  com¬ 
puter  simulation  to  generate  noisy  and  aberrated  tocai-piane  PSF  data.  A  now 
chart  ot  this  simulation  is  shown  in  Figure  3-i. 

CSDL  prepared  a  phase-retrieval  test  oased  on  aberrations  derived  trom 
cryogenic  deformations  ot  an  Itek  HALO  mirror.  This  mirror  is  circular  with  a 
diameter  ot  0.6  meter.  It  is  ultra-lightweight  and  made  trom  fused  silica. 

The  aberrated  PSF  was  sampled  by  an  8  x  8  array  ot  square  detector  elements, 
each  element  had  a  full  width  of  2.i3XF.  To  this  array  of  signals,  Gaussian 
random  noise  was  added  which  nad  a  unitorm  standard  deviation  ot  2  percent  ot 
the  peak  dif traction-limited  signal.  This  test  represented  Case  i. 

In  another  test,  Case  2,  PSF  was  decentered  such  that  the  peak  irradi- 
ance  of  the  dittraction-limited  spread  tunction  was  moved  by  (-0.984  XF , 

-0.984  XF)  which  corresponds  to  a  line-of-sight  error  of  (-0.984  X/D,  -0.984 
X/D),  where  D  is  the  diameter  of  the  system  pupil.  As  berore,  2-percent  noise 
was  added  to  the  sampled  data. 
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ABERRATIONS 


NOISY  DETECTOR  ARRAY  SIGNALS 


FFT  FAST  FOURIER  TRANSFORM 
PSF  POINT  SPREAD  FUNCTION 
OTF  OPTICAL  TRANSFER  FUNCTION 
DTF  -  DETECTOR  TRANSFER  FUNCTION 


Figure  1-1.  Flow  chart  of  computer  simulation  for  calculating 
noisv  and  aberrated  focal-plane  PSF  data. 
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In  the  third  test,  Case  3,  there  was  no  ime-ot-sight  error,  but  the 
noise  was  increased  to  3  percent.  The  three  8  x  8  arrays  of  aberrated  and 
noisy  signals  are  shown  m  Table  3-1. 

3 . 2 . 2 . 3  Preliminary  Test  Results 

By  using  their  algorithms,  Eikonix  estimated  the  aberrations  from  the 
detector  array  signals  in  terms  of  8,  i5,  and  23  Zernike  polynomial  coeffi¬ 
cients  that  did  not  agree  with  each  other  very  much.  They  also  differed  wild¬ 
ly  trom  case  to  case. 

To  compare  the  estimated  with  the  actual  aberrations,  CSDL  calculated 
aberrations  at  an  array  of  points  (within  a  circle)  from  the  Zernike  coeffi¬ 
cients.  When  the  estimated  aberrations  were  subtracted  trom  the  actual  ones, 
the  residual  aberrations  were  much  worse  than  the-  actual  ones  in  each  case. 
Generally,  the  standard  deviation  of  the  aberrations  increased  from  an  initial 
value  of  0.186  A  to  approximately  0.4  X. 

3.3  Image  Correction  by  Image  Sharpening 
3.3.1  Image  Sharpening  Technique 

The  image  can  be  corrected  m  a  closed-loop  manner  (as  opposed  to  the 
open-ioop  operation  of  the  phase-retrieval  technique)  by  optimizing  the 

sharpness  functions  obtained  rrom  the  focal-plane  data.  If  I(r)  represents 
the  focal-plane  image  distribution,  then  sharpness  functions 


S  ■  /  I^(r)  dr,  Extended  and  Point  Objects 

and 

«  /  I(r)  dr,  Poiat  Objects 

-► 

Ar 


attain  tneir  maximum  values  when  the  system  is  aberration  free.  The  first 
sharpness  function  uses  an  array  of  detectors,  but  the  second  uses  a  single 
detector  with  a  width  approximately  half  the  diameter  of  the  Airy  disc.  Small 
amounts  of  aberration  are  introduced  into  the  system  in  terms  of  Zernike  modes 
until  the  sharpness  function  is  maximized. 

3.3.2  Hardware  Demonstration 


The  image  sharpening  technique  was  demonstrated  using  an  adaptive  mem¬ 
brane  mirror.  Figure  3-2  is  a  schematic  of  the  laboratory  setup.  Six  Zernike 
modal  corrections  (detocus,  spherical,  two  astigmatisms,  and  two  comas)  were 
introduced  into  the  mirror.  Figure  3-3  is  an  example  of  image  correction. 
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Table  3-1.  Noisy  and  aberrated  detector  signals  sent  to  Eikonix. 

Case  I,  centered  PSF  and  2-percent  noise.  Case  2,  PSF 
decentered  by  (-0.9  84  XF,  -0.984  XF)  and  2-percent  noise. 
Case  3,  centered  PSF  and  5-percent  noise. 


Case  1 


8 

-18 

-9 

-5 

-1 1 

33 

13 

-7 

-15 

7 

6 

1 

7.1 

33 

30 

4 

-3 

0 

6 

-45 

-16 

2 

29 

179 

12 

27 

0 

5 

14 

-2 

6 

66 

(152} 

58 

-15 

-4 

4 

-2 

12 

1 

61 

95 

58 

-11 

3 

3 

29 

-20 

O 

16 

34 

i? 

5 

15 

2 

-2 

-8 

26 

35 

11 

-4 

-19 

-14 

1 

-16 

-8 

-12 

-8 

-15 

-30 

-6 

-25 

Case  2 

8 

26 

-3 

5 

44 

-13 

-4 

27 

-22 

7 

-22 

-33 

4 

21 

11 

5 

-8 

18 

6 

2 

13 

-11 

-29 

56 

14 

6 

9 

5 

-10 

3 

10 

102 

C240) 

29 

13 

1 

4 

-15 

-13 

41 

209 

164 

5 

-3 

21 

3 

-2 

11 

16 

-6 

31 

-2 

0 

16 

2 

-9 

31 

2 

-12 

20 

-3 

16 

-1 

1 

-22 

24 

-8 

-17 

-21 

-16 

-25 
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Case  3 

8 

-34 

-43 

14 

87 

-16 

-21 

51 

108 

7 

-30 

41 

22 

-9 

27 

22 

43 

-7 

6 

-5 

53 

-81 

-37 

166 

-27 

24 

-25 

5 

53 

-63 

7 

13 

29 

13 

40 

4 

57 

-2 

71 

15 

122 

-2 

14 

69 

3 

16 

2! 

-6 

-60 

56 

-41 

16 

-14 

2 

8 

-37 

12 

-2 

-1* 

-94 

-91 

-40 

i 

34 

-17 

-58 

45 

-49 

-41 

V 

-23 

-70 

Circular 

aperture 

center- 

to-center  spacing  and 

detector 

element 

si:e  *  2.13! 

in  XF-number  units.  Find 

the 

aberrations  in  terms  of  Zernike  coefficients. 

The  origin  (0,0) 

lies  at 

the  encircled 

number, 

and  the 

peak  diffraction- 

limited 

signal  is 

equal  to  835 

.  The  no 

ise  is  measured 

in  terras 

of  the  peak 

diffraction-limited  signal. 
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The  two  images  shown  In  this  tigure  represent  the  aberrated  and  corrected 
images  of  a  point  object. 
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Figure  3-2.  Schematic  ot  laboratory  setup  for  demonstration 
of  the  image  sharpening  technique. 


3.3.3  Software  Tests  and  Results 

Test  Case  1,  prepared  tor  Eikonix's  phase  retrip  al  test,  was  also  used 
on  the  CSDL  image-sharpening  algorithm.  Using  the  six  Zernike  modal  correc¬ 
tions,  the  Strehl  ratio  of  the  PSF  increased  from  0.26  to  0.62  in  one 
iteration. 

3.4  Review  of  HAC  Work  on  Image  Correction  Techniques 

Sam  Williams  and  his  associates  at  HAC  reviewed  their  work  on  OYSTER, 
color  algorithms,  image  momenta,  CORRWAVE,  and  phase  retrieval  at  DARPA  on 
13  Hay  1‘lfti.  The  phase-retrieval  work  was  presented  by  R.  Gonsalves  of 
Eikonix  as  a  subcontract  to  HAC. 

The  work  presented  consisted  of  computer  simulation  reeults  supported  by 
some  experimental  evidence.  They  felt  confident  that  these  methods  would  work 
under  appropriate  conditions.  Unfortunately,  these  conditions  were  not  des¬ 
cribed  adequately.  No  comparison  of  the  various  techniques  was  presented, 
every  technique  seemed  to  hold  promise.  How  a  particular  technique  would  work 
in  practice  was  not  discussed,  that  is,  system  related  issues  were  completely 
absent. 
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HAC  is  interested  m  developing  a  teat  bed  to  test  and  compare  their 
algorithms.  F*.om  what  has  been  reported,  it  seems  that  they  have  done  enough 
testing,  both  in  computer  simulations  and  in  the  lab,  to  show  that  these 
techniques  hold  promise.  These  techniques  should  be  reviewed  and  tested  by  an 
independent  party.  Per  ACOSS  li  Statement  ot  Work,  Paragraph  4. 1.4. 4,  CSDL 
will  prepare  test  cases  tor  the  HAC  color  algorithms  in  the  tr.ird  year  of  the 
three  year  program.  In  a  letter  to  RADC,  CSDL  recommended  the  following 
three-step  approach. 

■'  1 )  Prepare  and  Evaluate  Test  Cases  (Computer  Simulations)-  —To  carry 

out  this  task,  CSDL  can  work  with  HAC  in  the  same  way  as  work  has 
been  carried  out  with  Eikonix.  These  tests  will  bring  out  algo¬ 
rithm  subtleties  and  identify  limitations. 

(2)  Prepare  Deliverable  Algorithms — Depending  on  the  results  ot  Task 
1,  HAC  should  prepare  deliverable  algorithms  on  their  most  promis¬ 
ing  approach(es) . 

(3)  Algorithm  Hardware  Test — The  delivered  algorithms  should  De  tested 
in  a  hardware  simulation  such  as  OPTECAL.  Since  OPTECAL  is  en¬ 
visioned  as  a  system-level  optical-technology  simulation,  it  can 
test,  compare,  and  evaluate  component  technologies  and  algorithms.. 

3 . 5  Summary,  Conclusions,  and  Recommendations 

Table  3-2  summarizes  the  test  case  results  and  compares  the  phase- 
retrieval  and  image-sharpening  techniques.  It  is  evident  that  image  sharpen¬ 
ing  has  some  inherent  advantages. 


Although  Eikonix  is  still  working  on  the  te3t  cases  to  improve  the  algo¬ 
rithm  pertormance,  the  next  step  should  not  depend  upon  the  outcome.  As  a 
minimum,  a  hardware  simulation  should  be  prepared  to  test  the  practical  limits 
of  the  applicability  ot  their  algorithms.  CSDL  can  prepare  such  a  simulation 
using  their  in-house  image-sharpening  setup. 

The  adaptive  membrane  mirror  needs  to  be  repaired,  and  CSDL  has  the 
facilities  to  do  this  repair  and  make  the  mirror  operational.  Once  the  mirror 
is  operational,  the  hardware  facility  can  be  used  to  investigate  the  image- 
sharpening  technique  to  correct  aberrated  images  ot  extended  objects. 


Table  3-2.  Test  case  summary  and  comparison  of  phase-retrieval  and 
image-sharpening  techniques. 


Image-Plane 

Detector 

Object 

Operation 

Test  Case  Results 

Eikonix' s 
Phase 
Retrieval 

Array 

Point,  Yes 
Extended? 

Open  Loop 

Corrected  image  much 
worse,  a  increases 

from  U. Id  to  0.44  \ 

CSDL's 

Image 

Sharpening 

•  Single  for 
point  objects 

•  Array  for 
extended 
objects 

Point,  Yes 
Extended,  Yes 

Closed  Loop 

Six  Zernlke  modal 
corrections  increase 
Strehl  ratio  from 

0.2b  to  0.b2  with  one 
iteration 
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